Quantitative entanglement witnesses 
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Entanglement witnesses provide tools to detect entanglement in experimental situations without the need of 
having full tomographic knowledge about the state. If one estimates in an experiment an expectation value 
smaller than zero, one can directly infer that the state has been entangled, or specifically multi-partite entangled, 
in the first place. In this article, we emphasize that all these tests - based on the very same data - give rise to 
quantitative estimates in terms of entanglement measures: "If a test is strongly violated, one can also infer that 
the state was quantitatively very much entangled". We consider various measures of entanglement, including the 
negativity, the entanglement of formation, and the robustness of entanglement, in the bipartite and multipartite 
setting. As examples, we discuss several experiments in the context of quantum state preparation that have 
recently been performed. 

PACS numbers: 



INTRODUCTION 

Entanglement witnesses have proven tremendously help- 
ful in the experimental characterization of entanglement in 
composite quantum systems [1-13]. They are observables 
from the expectation values of which one can argue whether 
a prepared state is indeed entangled: whenever its expecta- 
tion value takes a value smaller than zero, then one can unam- 
biguously draw the conclusion that the state has been entan- 
gled in a particular fashion [1—4]: the entanglement has then 
been "witnessed". This approach seems particularly feasible 
or helpful in situations where one would like to avoid to col- 
lect sufficient data to arrive at full tomographic knowledge. 
Specifically in multi-partite settings when detecting multi- 
particle entanglement this can be costly. Also in instances one 
can tolerate larger errors when estimating entanglement wit- 
nesses compared to the procedure where one first estimates 
the full state. 

Originally, such a test for entanglement was thought to give 
rise to an answer to a "yes-no-question": the state is entangled 
or it is not. Yet, in this way, one does not make use of valuable 
information that one has collected anyway. Actually, one has 
implicitly recorded data that are sufficient to make a quanti- 
tative statement: if a test is very much violated - so delivers 
a value much smaller than zero - then one can infer that in 
quantitative terms, the state was highly entangled. This quan- 
titative statement is then meant in terms of some measure of 
entanglement. This is very useful information: One does not 
only know that the the specific entanglement property is con- 
tained in the state. But one can also give an answer to the 
question how useful a given state is, say, to perform a certain 
task of quantum information. 

This article emphasizes this fact, and advocates the 
paradigm of quantitative tests based on data from measuring 
witness operators Il4ll . Needless to say, one should under all 
circumstances only make use of the data that have in fact been 
acquired in an experiment, and avoid hidden assumptions con- 
cerning the nature of the involved states. But then, in turn, 



one should make use of the full information that can in fact be 
extracted from the measurement data, including quantitative 
assessments. 



PARADIGM OF QUANTITATIVE TESTS 

The paradigm we describe is the following: imagine one 
has collected data from a measurement of an entanglement 
witness, or a collection thereof. What is the worst case 
scenario one could have had, concerning the degree of 
entanglement? Certainly, one should provide conservative 
estimates in this context. This is typically the practically 
most relevant question: one has prepared a state, and wants 
to know to what degree one has succeeded in doing so. This 
test should make use of a minimal possible number of data, 
or measurement settings, certainly less than full tomography. 
So we aim for answers to 

"Given measurement data from measuring an 
entanglement witness, which one is quantitatively 
the least entangled state consistent with the 
data ? " 

This translates to an optimization problem: for a given mea- 
sure of bi-partite or multi-partite entanglement E, we aim at 
finding the solution of 



min E(p) 
subject to p consistent with the data, 



(1) 



or at least get reasonably good lower bounds. The general 
spirit of this paper will be to assume nothing more than the 
partial information provided by expectation values of entan- 
glement witnesses. Based on this information, we aim at find- 
ing good bounds to entanglement measures. We also comment 
on the tightness of these bounds. In fact, the provided strate- 
gies often give rise to the best (tightest) possible bounds based 
on this partial information. The "true state" of the system is 
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tx[Wp] = c 




FIG. 1: Schematic representation of state space. The set of separable 
states is depicted as the white region. The straight line represents 
an experimental test, so the hyperplane characterized by an entan- 
glement witness W and its expectation value tr[pW] = c. Then, 
one encounters a hierarchy of convex sets of states with increasing 
degree of entanglement, as quantified by any convex entanglement 
monotone. 



not assumed to be known, or it is not even assumed that it 
could in principle be measured, as full quantum tomography 
may be inaccessible. In turn, the optimization of entanglement 
witnesses, so the construction of tangent hyperplanes [11, 12, 
16-19] is an interesting (and computationally provably hard) 
problem in its own right, which we will not touch upon here. 
Any known findings in this field can however immediately be 
applied to our setting, in that the entanglement witness that is 
most violated will give rise to the best bound. We hence take 
the entanglement witness as such and the corresponding data 
for granted, and will provide good bounds for entanglement 
measures based on them. This is actually the situation one 
faces when interpreting experimental data. 

There is a body of work somewhat similar in spirit in the 
literature. The need for conservative estimates, so for mini- 
mizing the degree of entanglement in the context of a Jaynes 
statistical inference scheme consistent with the data was al- 
ready noted in the early work Ref. B20I1 . Also, conceptually, 
this is related to a connection of violations of Bell inequalities 
to entanglement measures 12111 . and to quantum state estima- 
tion as in Ref. [22 ] . 

In this work we consider the bi-partite and multi-partite 
setting. The system can hence be thought of consisting of a 
number of subsystems, such that the Hilbert space is given by 
H = € dl ® ■ ■ ■ ® € dN . We assume that we have collected 
data that we can estimate {ci, . . . , c„} based on a number of 
entanglement witnesses Wi, . . . , W n , meaning that 



(Wi) = tr[Wi/>] 



(2) 



for i = 1, . . . , n, for example, for a single entanglement wit- 
ness, n = 1. Witness W means in the bi-partite setting, 
N = 2, that for all separable states 



P = 



E 



pip 



(i) 



p 



(2) 



(3) 



on H = € dl ® € d2 we have that [1-4] 
tr[Wp] > 0, 



and at least for a single entangled state p, one finds that 
ti[Wp] < 0. 



(5) 



This is very intuitive: the separable states form a convex set, 
and the witness defines a hyperplane in state space that sep- 
arates the separable states, see Fig. Q] In the same way, one 
can define entanglement witnesses for the various classes of 
multi-particle entanglement, in a setting with Hilbert space 
H = <D dl <%) ■ ■ ■ <%) (D djv . For witnesses in infinite-dimensional 
systems and relationships to entanglement measures, see Refs. 
[23-26]. In this paper, we refrain from introducing these 
multi-partite entanglement classes, and refer for that to Refs. 
[27-30]. 

We should mention at this point that if one allows for wit- 
nesses taking several identically prepared specimens into ac- 
count, one can often improve the bounds to entanglement 
measures. On the positive side, this gives rise to sharper or 
tight bounds, often making use of few different types of mea- 
surements, or indeed even single ones [31-35]. On the neg- 
ative side, one needs to implement collective operations, ei- 
ther with quantum networks, or in optical settings, with joint 
operations involving bringing together independent sources 
at beam splitters. Nonetheless, the first experimental mea- 
surements of a two-copy witness for arbitrary two qubit pure 
states was recently reported 13511 . Although we make the pre- 
sented ideas explicit for the most frequently applied approach 
of measurements on individual specimens, it should be noted 
that many of the presented ideas are also applicable to this 
case of collective operations. 



IMPLICATIONS TO ENTANGLEMENT MEASURES 

Subsequently, we will present a framework of quantitative 
tests, and discuss a number of bounds for different entangle- 
ment measures. We will also discuss several examples, taken 
from the bi-partite and multi-partite context. The concept of 
the conjugate function will play a central role here. 



Negativity 

In this first part we consider bi-partite splits of our sys- 
tem: so the system is either naturally bi-partite, or we group 
the subsystems into two parts, with joint Hilbert space H = 
<C dl (g>(D d2 , with state space S(H.). The negativity is a measure 
of entanglement defined as 



En(p) 



\p% 



1. 



(6) 



in terms of the trace-norm ||j4||i = tr|A|. p v denotes the par- 
tial transpose of p. The negativity has been introduced in Ref. 
03], compare also Ref. Q7h . and independently shown to be 
an entanglement monotone in Refs. 1 38, 3^1 . The logarithmic 



(4) version log 2 ||p r ||i of the negativity is also an entanglement 
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monotone |40 L and a useful upper bound to the distillable en- 
tanglement |39l 41]. For this measure of entanglement, we 
will indeed find very simple, yet tight and useful bounds. 

What we are interested in here is the minimally entangled 
state consistent with what has been measured. So we seek the 
solution of 

E NMn = in{ H/lli-1 (7) 
subject to tr[pWi] = a, 

tr[p] = 1, 

which is the desired quantity. Now, this can also be written as 

En, min = inf maxtr[Pp r ] - 1 (8) 
subject to ||-P||oo = 1) 
tr[pWi] = a, 
P>0, 
tr[p] = 1, 

as \\A\\i — maxtr[AA] with a maximation over all operators 
with 1 1 -X" ||oc = 1, according to the variational characterisa- 
tion of the trace-norm. In turn, obviously, any such A with 
1 1 X I |oo = 1 gives rise to the lower bound 



-EiV.min > inf 

subject to 



tr[Ap r ] - 1 
tr[pWi] = Ci, 
P>0, 
tr[p] = 1. 



(9) 



We can now take an operator consisting only of the partial 
transposes of the witnesses we have measured, 



A' 



am 



1. 



Oj 6 E for i = 1, . . . ,n + 1, such that ||X||oo = 
there is nothing to minimize any more, as tr[p r A r ] 
we arrive at 



> 2J ctid + a n+1 - 1. 



(10) 

1. Then, 

= ti[pX]: 

(11) 



This is indeed a very simple bound. Yet, it is a useful, and 
tight one. 

How can one find a suitable choice for ai, . . . , a n +i? Any 
choice such that — 1 < X < 1 as in Eq. ( fTot gives rise to 
a bound. In turn, one can also find the optimal choice in an 
efficient manner: The problem we encounter is, 



subject to 



^ aid + a n +i 
i=i 

-t<X <1, 

n 

A = 



(12) 



C*ri+ll, 



as an optimization problem over {a±, . . . , a n +i}- This is an 
optimization problem we can run beforehand: it is actually 
a semi-definite optimization problem 14211 . so an optimization 
problem that can be efficiently solved, with certifiable error 
bounds. But for the use of this criterion as such, one does no 
longer have to solve any optimization problem. 

At this point, a remark is in order concerning the tightness 
of the constructed bounds. Let p T be the partial transpose of a 
state on <C dl 55 <C d2 , N p eigenvalues of which are strictly posi- 
tive, N n eigenvalues are strictly negative, and Nq eigenvalues 
take the value 0. Then any A satisfying 



l/lli 



tr[Ap r 



(13) 



has a spectrum containing at least N p times the value 1 and 
N n times the value —1. In our present context, this means that 
for a given system dimension <C dl ® (D d2 , the above bound 
Eq. dTTb is tight whenever the there exist states p such that 
S p > N p and S n > N n , where S p and S n are the number 
of ±1 eigenvalues of A, respectively. Then the bound is just 
saturated by actual physical states. 

Example 1 (Bound to the negativity) As a very simple ex- 
ample, consider states on (D 2 ® (D 2 . The witness we take is 



Wi = \<r)(<t> 



(14) 



which is an optimal entanglement witness, in that it is tangent 
to the set of separable states. Here and in the following, |^> ± ) 
and denotes the state vectors of the familiar Bell states 
for two qubits. Now consider A = —2W[ + 1, so a\ = —2 
and«2 = I- The matrix A clearly satisfies ||A||oo = 1- Then, 
whenever we get a value tr[Wip] = c, we can assert that 



E N (p) > 2\c\. 



(15) 



It is also easy to see that this bound is tight: The spectrum of 
A is given by {1, 1, 1, —1}. A family of states saturating the 
bound is given by 



Ai^ + )(v; + i + (i-A)i<n(vn 



(16) 



for which E N {p) = 2\c\ = |2A - 1|. For c = -1/2, the 
only state consistent with this value is the maximally entan- 
gled state \4> + )(i>+\, yielding E N (\tp+)(ip+\) = 1. 

In turn, we can see what we may gain from using two wit- 
nesses: 

Example 2 (Bound from two witnesses) Let us take the two 
entanglement witnesses Wi — \<fii) (4>i\ T , i = 1,2, where 

l^ = >,o) + l|o,i) + i | i ) o) + (|) 1/2 |i,i) ) 
Ifc) - >.o> + ^lo.i> + i|i,o> + (^) 1/a |i,i>, 

(17) 
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and ci = —1/3 and c 2 = —1/6. Then, we may evaluate 
the optimal bound based on each witness separately, and the 
best bound based on both simultaneously. From solving the 
semi-definite optimization problem, we find in case of W\, 



E 



./V,min 



then for Wi, 



E 



iV.min 



>2/3, 



> 1/3. 



(18) 



(19) 



Indeed, in the combined case using W\ and W2, we obtain the 
better bound 



E NMn > 0.7375. 



(20) 



This shows that the suitable processing of several witnesses at 
the same time can give rise to optimized bounds. The bound 
arising from the data from two witnesses is stronger than each 
bound resulting from either of them. 

The presented bounds are based on simple witnesses for 
qubit systems, but it should be clear that the construction is 
general enough such that bounds can be identified in fact for 
arbitrary entanglement witnesses in any dimension. 



again S(H) denoting state space, which is 

f*(X) = sup {tr[|^)(V|X]-/(|V)(V|)} (26) 
Wen 

for concave functions /. In turn, the conjugate function of the 
conjugate is the convex hull of the function itself 14211 : In other 
words, since the entanglement of formation is the convex hull 
of the reduced entropy function itself, we have that 



where 



rip) = hp), 



r(p) = S up{tr[px}-r(x)}. 



(27) 



(28) 



By definition, f*(X + at) = f*(X) + a for any X. Now, 
for any 



X = Y J OL i W i , 



(29) 



i=l 



we indeed arrive at the bound 



Convex hull measures and the conjugate function 

Many entanglement measures are defined as a convex hull 
of a function, so as / = cof. This is nothing but 



hp) = mm \ ^Ptfipi) ■■ Yl 

\ i i 



PiPi = p 



(21) 



for states p. The most familiar example of this sort is the 
entanglement of formation, for which this function / is the 
reduced entropy function 



/(p) = (Sotr 2 )(p), 



(22) 



where tr2 is the partial trace in a bi-partite system and S(p) = 
— tr[plog 2 p] is the von-Neumann entropy. The convex hull of 
a function / can alternatively also be written in the form 



t{p) = sup{tr[Xp] : V|V) € H : 

){m < fi\i>){ip\)}. 



X 

till 



(23) 



Note that we have for consistency assigned f(x) = 00 in case 
of x < 0. Again, we aim for bounding the solution of 

S min =inf /(p), (24) 
subject to tr[pWi] = Cj, 
P>0, 
tr[p] = 1, 

i = 1, . . . , n. We can make use of the conjugate function [42- 
44], also known as the Legendre transform: This is defined as 

fix) = sup {tr[pX]-/(p)}, (25) 



(30) 



in terms of the conjugate function /* of /. In this way, we do 
not have to evaluate the convex hull explicitly. 

Moreover, the bounds constructed in this way are always 
tight. It follows from the duality of the convex hull of the 
function and its Legendre transform that the bounds are tight 
when varying over all {ai, . . . , a„}. There always exists a 
state p satisfying E m [ n = /**(p), so for example Ejr^m = 
Ep(p) for the entanglement of formation. In this sense, the 
given bounds are the best possible bounds of this form. 

In case a symmetry can be identified, the estimation of the 
conjugate function of a given function can be simplified. To 
bring the conjugate function into a form that is more accessi- 
ble to numerical assessments, we can proceed as follows: If 
/ = id tr 2)> and g is concave, we can define 

g(p) - MMYp]-g'(Y)}, (31) 

g'(Y) = in£MYp}-gip)} > (32) 
p 

and can write the above conjugate function as (assuming that 
X and Y\ are Hermitian) 

fiX) = supsu P {tr[(X-(y 1 ®l))p]+ 5 '(r 1 )} 

P Yi 

= sup{A max (X-(Fi®l)) + g / (yi)}. (33) 
n 

For the entropy function g(x) = —x log 2 x < f° r example, the 
conjugate g' is known, and one finds 142 ] 



5 '(y 1 ) = -log 2 tr[exp(-F 1 log2)]. 



(34) 
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In this form, the problem is in a suitable form for such numer- 
ical assessments. The resulting bound is then a combination 
of the numerically evaluated expression and the value for c 
from the actual data. In practice, this numerical evaluation 
amounts to a global optimization problem, which can, for a 
small number of parameters in typical problems in the quan- 
tum information context, be solved for an arbitrary witness. 
Also, semi-definite relaxations as in Refs. [ IT , 46| readily give 
rise to certifiable bounds. 

As an example, let us look at the entanglement of formation, 
and a single witness W\ . Then, 



-E^min > ol\c- f(aiWi), 



(35) 



with a,\ G R, so any choice for u\ delivers a bound. Obvi- 
ously, an optimal bound is achieved using 



E F ,mm > sup ai {ajc - f*{a 1 W 1 )} 



(36) 



Similarly, more than a single witness can be considered. So 
one needs to find good upper bounds to f*(aiWi). 

Example 3 (Bound to the entanglement of formation) 

This becomes particularly simple for witnesses of the form 
Wi = |0)(0| r in C d <g> € d , for entangled state vectors \4>). 
We consider the conjugate function /*, evaluated at a{W\. 
It is easily seen from Eq. (23) that, for any entanglement 
measure, the conjugate function is invariant under local 
unitaries. Then, without loss of generality \<j>) can be taken to 
be of Schmidt form 



d 

E 

i=i 



&|m)- 



(37) 



The partial transpose |</>)(^>| r gives rise to the form 

d 

\4>)(4>\ T = E O&b'^XM, (38) 

so in a product basis a direct sum of 1 x 1 and 2x2 matrices. 
We seek the maximal value of aitrfW^i (ip\] and a minimal 
value for / ( | ip) (tp | ) . Let 



■ , nlax .. - 



(39) 



sow = £k£,i for some k, I = 1, . . . , d. It is not difficult to see 
that then the optimal state vector takes the form 

= a\k,l) - (1 -a 2 ) 1/2 |Z,fc) (40) 

for some a <G [0, 1]. This state vector gives rise to 

tr[-Wi|V>)(^|] = 2a(l -a 2 ) 1/2 w (41) 



and (^|) 



-a 2 log 2 a 2 - (1 - a 2 )log 2 (l 



defining the concave (classical entropy) function h{x) 



— x\og 2 x — (1 — x )l°g2(l — x ) anc l tne concave function 

g(x) = (x{l - x))V\ 



f*( ai Wx) = sup <^ -2a!wa(l - a 2 ) 1 ' 2 

oG[0,l] I 

+ a 2 log 2 a 2 + (1 - a 2 ) log 2 (l - a 2 



= sup {-2a 1 wg(p) - h(p)}. (42) 

P e[o,i] 

We can distinguish three regimes. Define the parameter b = 
—2a\w and the function z(p) = bg(p) — h(p). The second 
derivative of z is given by 



z"(p) 



1 



log(2)p(l-p) 4(p(l-p))3/ 2 " 



(43) 



The function z is convex iff z" is non-negative for all p £ 
[0,1], which occurs when 



4(p(l -p)) 1/2 

b < min — = 0. 

" P log 2 



(44) 



The function is concave when z" is non-positive for all p, 
which occurs when 

b > max4p(l -p) 1/2 /log(2) = 2/log(2) = 2.88539. 
v 

(45) 

In between these values, z is neither convex nor concave. If z 
is convex (b < 0), its supremum occurs at one of the extreme 
points, either p = or p = 1. But of course, either one gives 
the same value, namely 0. If z is concave (b > 2/log(2)), 
it has one supremum. By the even symmetry of z around 
p = 1/2, the supremum must occur at p = 1/2, yielding 
as supremum value (b — 2)/2. For determining the supremum 
in the case < b < 2/ log(2), a transcendental equation has 
to be solved. The supremum as function of b can be approxi- 
mated from above by the polynomial 



0.0018766 + 0.00823% 2 + 0.0197336 3 
- 0.0056496 4 + 0.0014306 5 . 



(46) 



The average error of this approximation is 0.00017. We may 
take ai = —1/w, then /*(aiWi) < c with c = 0.14985. 
Therefore, we obtain the bound 



E 



.F,min 



> 



\c\/\w\ - cq. 



(47) 



To emphasize that again, we do not assume the "true state" 
to be detected to be known or accessible. For completeness, 
we do elaborate on an example showing the tightness of the 
bound. For example, for the family of states in Eq. ( TToT ) for 
A G [1/2, 1], and for the witness W\ = \<j>~) (cf>~ | r as in Eq. 
CLUl, we find 2\c\ = |2A- 1|, w = 1/2, and in fact 



E F (p) = h(l/2+(\(l-\)) 1 / 2 ) 



(48) 
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which has to be compared with 



for q > 1, we get an upper bound to /* (X) as 



E F , mm > |2A- 1| -co- (49) r(X)<sup{tr[|^)(^|X]- 5g (|^)(V|)}. (57) 
As can easily be seen, this is a very good lower bound (as a 

tangent the best possible affine bound), and the bound is tight The function 9i is no longer concave, but we nevertheless get 

for A = 7056 an a PP ro P r i ate bound when optimizing over pure states. Par- 
ticularly useful is the case of q = 2, when we merely need to 

Example 4 (Second bound to the entanglement of formation) evaluate tr [tr 2 [ | ip) ( -0 1 ] 2 ] ■ 
Let us consider a witness of the common form 



Wi = at-b\(f>}((f>\ 



(50) 



with a, b > 0. The conjugate function can be easily written as 

f*(a l W 1 ) = sup {aio-oi6|^|V)| a -/(|V>W)}-(51) 
men 

We assume \<f>) to be in its Schmidt form, given by Eq. d3Tb . 
The state vector | ip) might also be written as 



(52) 



where the basis {\i')} is not necessarily equal to the Schmidt 
basis {\i)} of \<p). One can thus minimize over the basis {\i')} 
and the Schmidt coefficients /i;. The last term in the right hand 
side of Eq. ( BTT l clearly does not depend on In turn, 

given a fixed set of Schmidt coefficients fa, this implies that 
the optimal basis will be the one which maximizes the overlap 
| ((f> | if>) | 2 ■ It can be easily shown that the maximum is obtained 
when choosing {\i')} to be equal to the Schmidt basis of \<f>). 
Therefore, we are left with an easier maximization problem, 
over the Schmidt coefficients only, given by 

f*(a 1 W 1 ) = supjaia- aib^Y^&fiij (53) 
+ & 2 log 2 (M 2 ) 

i 

Although it is not possible to solve Eq. d53l analytically in 
terms of ai,a and b for all choices, it can be easily numeri- 
cally evaluated. For example, let us consider a = a± = 1, 
6 = 3/2, and 



1/2 



10,0} 



1/2 



IM>. 



(54) 



so & = (1/3) 1 / 2 and £ 2 = (2/3) 1 / 2 . We then get the follow- 
ing bound for the entanglement of formation 



E 



F'.min 



> Id 



0.5550. 



(55) 



Example 5 (Bounds from Renyi entropies) Since for the 
entropy function /, we have that 



f{x) > g q (x) 



1 



iog 2 (^) 



(56) 



Further bounds can, e.g., be found in Refs. 1471 14811 



Remarks on exploiting symmetry 

If one has a witness which is invariant under a local sym- 
metry group, one can in instances simplify the evaluation of 
good bounds under the constraint provided by the entangle- 
ment witness: One can take the Haar average with respect to 
that group, which will always diminish the degree of entan- 
glement. So a twirling with respect to a, for example, U (£> U, 
U <g> U*, or O <£> O-symmetry, or one corresponding to SU (2) 
or symmetric group representations, can only give a lower 
bound [49]: For any convex entanglement monotone /, 



E min = inf 
subject to 



f(p) 

tt\pWx] = ci, 

p is a symmetric state. 



(58) 



Hence, we have to evaluate an entanglement measure under 
symmetry [49-53], given the constraint. 

Example 6 (Symmetry) To give a very simple example, let 
us consider a witness of the form 



W x =al + b\(f>-)((f-\ l 



(59) 



for some a, b g R and for states on (D 2 <£> (D 2 . Since \<f>~) (<f>~ \ 
is a J7 <£> CZ-symmetric state, this witness is U <E> J7*-symmetric. 
Therefore, we can optimize the bound with respect to U ®U*- 
symmetric states, which is the one-dimensional convex set 



P = Al/4+(l-A)|^+)(V; + |, 



(60) 



for A G [0, 1]. The entanglement of formation Ep(p) of such 
symmetric states p, in turn, is known 114911 . 



Concurrence 

An interesting example where the conjugate function can be 
analytically calculated is the concurrence of two qubits 154 1. 
Let us define the following basis for (D 2 (g> <D 2 , 



|*o) = 77|(l > > + I 1 . !»» 1*0 = -7f(l°'°) - I 1 ' 1 »' (61) 
|* 2 ) = -J=(|0, 1) + |1, 0)), |* 3 ) = -^(|0, 1) - |1,0))(62) 
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As is well-known, a general two-qubit pure state can then be 
written as 



I0W> = 

The concurrence of a pure state is defined as 



3 

£■ 

i=0 



3 
i=0 



(63) 



(64) 



and extended to mixed states by a convex hull construction. 
The importance of the concurrence is twofold. On one hand, 
it is intimately related to the entanglement of formation of two 
qubits. Indeed, given the concurrence of a two qubit state p, 
its entanglement of formation reads 15411 



E F (p)=H(± (l + (l-C(p) 2 ) 1/2 )) 



(65) 



On the other hand, an analytical expression for the concur- 
rence of a general two qubits mixed state is known 1 54], which 
in turn implies an analytical formula for the entanglement of 
formation. 

The first interesting bound to the concurrence based on a 
witness, derived in Ref. 115 511 . is given by 



C(p) — max < 0, 



min tr[M) (AT pi 

AeSL(2,C) 



(66) 



where \A) denotes the unnormalized state vector (A <£> 1)|J) 
with 1 7} = |0, 0) + 1 1,1) and A is any (in its determinant 
normalized) 2x2 invertible matrix. It is thus seen that any 
witness Wi of the form 



W 1 = \A)(A\ l 



(67) 



provides a lower bound to the concurrence. Note that this class 
of witnesses is exactly the class of optimal entanglement wit- 
nesses of two qubits, where optimality refers to the robustness 
with respect to white noise (mixing with the identity). 

Although Eq. (l66l l constitutes a useful tool to estimate C, 
it has the drawback that only witnesses of the restricted form 
given by Eq. d67| i can be used. On the other hand, the method 
based on conjugate functions, oulined in the previous subsec- 
tion, can be applied to any entanglement witness. We now aim 
at showing that the conjugate function of the concurrence can 
be evaluated analytically. 

The concurrence can be expressed as the convex-hull of the 
function 

/(p) = (2(l-tr[tr 2 [p] 2 ])) 1/2 , (68) 

defined on states p. It is easy to check that / is concave. This 
in turn implies that the supremum in Eq. ( f26b can be calcu- 
lated over pure states only. The conjugate function /* can be 
expressed as the following optimization problem, 



f*(X) 



subject to 



Ci\^i) is normalized. 



The optimal solution, as a function of X, although this not 
being a convex problem, can be readily evaluated with the help 
of a computer algebra program [56]. 



Further convex roof measures 

Note that we considered the entanglement of formation as 
an example for a "convex roof measure". There are other im- 
portant measures of entanglement in the multi-partite context, 
where the presented ideas can be applied. This applies in par- 
ticular to the geometric measure of entanglement. This is a 
measure for entanglement in the multi-partite case, which is 
defined for pure states as 115 711 



inf 

p 



(70) 



where ||^4||2 = tr[v4 2 ] is the Hilbert-Schmidt norm, and the 
infimum is taken with respect to all pure product states. The 
extension to mixed states is done via a convex roof construc- 



tion. Similarly, the global entanglement of Ref. 15811 may be 
considered. Both quantities are proper multi-partite entangle- 
ment monotones. For a survey on multi-partite entanglement 
measures, see, e.g., Refs. 129115911 . 



Robustness 

Given a bi-partite or multi-partite state p, its generalized 
robustness of entanglement - introduced in Refs. |6(J 61] - is 
defined as the minimal s > such that the state 



(71) 



l + .s 



is separable, where a is another arbitrary state. This measure 
can be interpreted as the minimum amount of noise neces- 
sary to wash out completely the quantum correlations initially 
present in the state p. In addition, the generalized robustness 
also has the operational interpretation for bi-partite systems 
as the usefulness of the state in question as an ancilla in tele- 
portation protocols 116211 and is a multi-partite entanglement 
monotone loHl . 

For our purposes, a very convenient representation of the 
gen eralized robustness, obtained as the Lagrange dual form 
[il of Eq. (EB |63|. is 



Er{ P ) 



: jo, - rmntr[Wp] j 



(72) 



where W is varied over the set of witnesses with maximum 
eigenvalue smaller than unity (W < 1). Note that by con- 
sidering different sets of witnesses, one can quantify all the 
different kinds of multi -partite entanglement [27-29]. 

As discussed in Refs. 163l 16411 . it follows directly from Eq. 
(1731 that the expectation value of any measured witness W 
gives rise to a useful lower bound to the generalized robust- 
ness. Then, when tr[Wp] = c < 0, 



E RMa {p) > \c\/X max (W). 



(73) 
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In full generality, this same approach can be applied to any 
entanglement measure which can be expressed as 



E(p) = max < 0, — min tr[Wp] 



(74) 



where M is the intersection of the sets of entanglement wit- 
nesses with some other set (e.g., the set W < 1). Inter- 
estingly, several other well-known entanglement quantifiers, 
such as the best separable approximation 16511 . the Rains fi- 
delity of teleportation l66Tl . and the concurrence (see Eq. 
d64ll), fit into this classification. Here for concreteness we 
focus on the generalized robustness and on the random ro- 
bustness, which we discuss in the sequel. 

A source of noise often considered in experiments is so- 
called white noise, in which the initial state p is driven to a 
state of the form 



1 

P + S D> 



(75) 



where s is related to the amount of noise introduced in the 
system. Here D stands for the dimension of the Hilbert space 
which p acts on. In this sense, it is interesting to ask what 
is the maximal tolerance of an entangled state to white noise, 
before all its initially entanglement is transformed into merely 
classical correlations. The random robustness l60ll is exactly 
such a quantity. In the framework of Eq. d74b , we can express 
it as the minimization over the set of entanglement witnesses 
with trace equal to D. Hence, every entanglement witness W 
can be used to lower bound it as 



E r ,min(p) > 



D\c\ 

ti[wy 



(76) 



again with tr^p] = c. 



Example 7 (Tri- and quadripartite photonic entanglement) 

As an example, we consider two multi-partite witnesses 
which have been measured in the photonic parametric-down- 
conversion experiment of Ref. [5], so A = 3 and A = 4. 
Consider the following multi-partite pure states vectors 

\W) = -^=(|0 5 0, 1) + 10,1,0) + |1,0,0» (77) 



and 



|^ 4 >H-±=(|0,0,1,1) + |1,1,0,0) 



(78) 



- -(|0, 1, 1,0) 4- 11, 0,0, 1) + |0, 1,0, 1) + |1, 0, 1, 0))). 

Then the two associated multi-partite entanglement witnesses 
which have been measured are given by 



W w = ^l-\W){Wl 

w vw = jji-|* (4) )(* (4) | 



(79) 
(80) 



Whereas the witness Ww detects genuine tri-partite entangle- 
ment, having positive values on separable and bi-partite entan- 
gled states, the operator W^,(4) witnesses genuine four-partite 
entanglement, being positive on separable, bi-separable, and 
tri-separable states. The measured expectation values in turn 
are J5J], 

tx[W w p] = -0.197 ± 0.018, (81) 
tr[W<pC4)/o] = -0.151 ±0.010. (82) 

Hence we readily have the following estimates on the robust- 
ness of the, a priori unknown, measured states p\ and p2, con- 
siting of three and four parties respectively: 

E RMn ( Pl ) > 0.2955 ±0.027, E rM n(pi) > 0.360 ± 0.096, 
E R , mia (p 2 ) > 0.201 ±0.013, E^ mn {p 2 ) > 0.220 ±0.021. 

(83) 

Example 8 (Four-photon graph state) In Ref. |(J, entangle- 
ment witnesses have been employed to characterize optical 
four-photon graph states [70-74] that have been prepared from 
entangled photon pairs, followed by a controlled-phase gate 
(compare also Ref. [69]). For the four-photon cluster state 
1 7511 . N = 4, the given witness is 



W CW = 31 - i (Z<»ZV> + l) ( Z< W 4 > + 1 



- 1 (XW»ZW + l) (Z^Z^ + l) . (84) 

The maximal theoretical value is tr[W c <4) p] = — 1, the mea- 
sured value is 



tr[W c( 4)p] = -0.299 ±0.050. 
This gives rise to 



E 



> 0.0997 ±0.0167, 
> 0.1120 ±0.020. 



(85) 



(86) 
(87) 



Example 9 (Quantum byte) In the recent spectacular exper- 
iment of Ref. |7J], 8 ions have been prepared in a multi -particle 
entangled state. The multi-particle entanglement has in turn 
been demonstrated using the concept of entanglement wit- 
nesses. In order to introduce the multi-partite entanglement 
witnesses that have been measured, we have to consider the 
A-partite W states 

\W N ) = (|0,...,0,0,1) + |0,... ,0,1,0) (88) 
+ |0,...,1,0,0) + --- + |1,...,0,0,0))/VA. 

Define the A-qubit state vectors \BSi) = |A) ® \Wn-i), 
which consist of |0) on the i-th qubit and the state vector 
|Wjv-i) on the remaining qubits, and the corresponding op- 
erators 

JV 

Q N = 10\Wn)(W n \ ~PnJ2 \BSi){BSi\, (89) 
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where /3jv is a fixed real number for each value. Next, define 

7tv = max|^) = | )g)| 6 )(*|Q|*), where ranges over all 



possible bi-separable state vectors II27L 12911 with respect to all 
possible bi-partitions. The witnesses are then given by 



Wi 



Byte 



7ivl- Q 



N- 



(90) 



They hence classify tri-partite entanglement. As explained in 
Ref. 10], the expectation values reported refer to the normal- 
ized versions of Weyte with trfWsyte] = 2^, where N is the 
number of parties of the state. Therefore, we can readily read 
them as lower bounds to the random robustness of the un- 
known 3, 4, 5, 6, 7, and 8 qubit states: 



E r 
E r 

E r 



,(3) > 0.532, £ r . mm (4) > 0.460, (91) 
,(5) > 0.202, E rMa (6) > 0.271, (92) 
,(7) > 0.071, E rMn {8) > 0.029. (93) 



Needless to say, the same discussion can also be performed 
based on other multi-partite measures of entanglement, such 
as the geometric measure of entanglement. 



CONCLUSIONS 

In this work, we have introduced quantitative bounds to en- 
tanglement measures, on the basis of expectation values of 
entanglement witnesses. In this way, quantities that are fre- 
quently measured in order to detect entanglement in experi- 
mental settings can be augmented with a stronger, quantitative 
statement on the degree of entanglement. In most instances, 
this does not require any additional effort at all, but these 
quantitative bounds may even be added in retrodiction. Sev- 
eral measures of entanglement have been considered. Need- 
less to say, similar methods can also be made use of if one 
has additional knowledge at hand about the system, say, from 
correlation measurements. 

We have considered the concept of the conjugate function 
in this context, and have presented a number of new bounds 
to entanglement measures. For the negativity, simple and very 
useful bounds emerged from a variational principle. We have 
discussed several examples taken from experimental settings, 
both from the context of linear optics, as well as of trapped 
ions. In this way, we have sharpened the notion that if a test 
for entanglement is indeed violated to a large extent, then the 
degree of entanglement can be expected to be large. 
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